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1. 


INTRODUCTION: 


High  grade  serous  ovarian  carcinoma  is  among  the  most  lethal  cancers  affecting  women  in  the 
U.S.  While  most  therapeutic  approaches  have  focused  on  malignant  epithelial  tumor  cells  and 
their  genetic  alterations,  it  is  becoming  increasingly  clear  that  the  tumor  microenvironment  plays 
an  equally  important  role  in  tumor  evolution.  The  presence  of  cancer  cells  induces  a  reaction  in 
the  surrounding  stromal  cells  similar  to  fibrosis  after  an  injury.  These  reactions  can  also  reduce 
therapeutic  efficacy  of  chemotherapy  by  creating  a  physical  barrier  for  drug  transport  while 
providing  a  protective  environment  for  cancer  cells  to  repopulate  after  completion  of  treatment. 
Thus,  it  is  thought  that  anti-cancer  therapies  should  target  not  only  malignant  cancer  cells  but 
also  the  microenvironment  that  fosters  tumor  growth  and  survival.  Our  goal  is  to  demonstrate 
that  targeting  processes  responsible  for  the  formation  of  carcinoma  associated  fibroblasts  (also 
known  as  CAFs)  in  the  tumor  microenvironment  will  effectively  attenuate  tumor  growth, 
improve  intratumoral  drug  delivery  and  restore  anti-tumor  immune  responses.  We  are  using  three 
different  approaches  to  targeting  CAFs.  The  first  approach  is  to  We  will  test  anti-fibrotic  agents 
for  their  efficacy  in  preventing  CAF  activation  and  increasing  sensitivity  to  chemotherapy  in  a 
mouse  model  of  ovarian  cancer  that  was  developed  in  our  laboratory.  The  second  approach  is  to 
increase  the  precision  of  targeting  activated  CAFs,  by  targeting  a  protein  that  we  previously 
identified  to  be  present  in  activated  CAFs  but  absent  from  fibroblasts  associated  with  non- 
cancerous  conditions  such  as  fibrosis,  inflammation,  and  wound  healing.  The  third  approach  is 
test  several  agents  for  their  efficacy  in  inducing  CAF-to-cartilage  differentiation  with  the  idea 
that  a  terminally-differentiated  microenvironment  cannot  protect  malignant  cells  from 
chemotherapy  or  foster  their  dormancy  for  future  recurrence. 

2.  KEYWORDS: 

Ovarian  cancer,  tumor  microenvironment,  cancer- associated  fibroblasts,  fibrosis,  targeted 
therapy,  clinical  outcome 

3.  ACCOMPLISHMENTS: 


What  were  the  major  goals  of  the  project? 


Specific  Aim  1  (specified  in  proposal) 

Timeline 

Percent 

Completed 

Major  Task  1 

Test  the  therapeutic  efficacy  of  CTGF,  CTSK,  FN1,  and  LOXL2 
inhibitors 

Months 

Cedars-Sinai 

Medical 

Center 

Subtask  1 

Amend  approved  IACUC  protocol  5318  (Mouse  Models  of  Tumor 
Microenvironment,  PI:  Orsulic)  for  local  approval  and  send  related 
material  for  DoD’s  approval. 

Upon  award 
notice 

100% 

Subtask  2 

Purchase  FVB  mice,  drugs,  and  reagents;  plan  experiments. 

1-2 

70% 

Subtask  3 

Implant  FVB  mice  with  mouse  ovarian  cancer  cells. 

2-25 

25% 

1 


Subtask  4 

Treat  mice  with  CTGF,  CTSK,  FN1,  and  L0XL2  inhibitors 

Assess  therapeutic  efficacy: 

1.  Tumor  growth:  tumor  weight/volume,  luciferase  whole-animal 
imaging. 

2.  Tumor  invasion  and  metastasis:  dissection  and 

immunohistochemistry. 

3.  Stromal  differentiation:  Masson’s  trichrome  stain,  qPCR  and 
immunostaining  for  myofibroblast  markers  (a-SMA,  fibronectin, 
COL11A1). 

4.  Chemotherapy  diffusion:  quantification  of  fluorescently-labeled 
dextran  beads. 

5.  Tumor- infiltrating  immune  cells:  flow  cytometric  analyses  with 
antibody  cocktail  (CD3,  CD4,  CD8a,  CD44,  CD62L,  CD25, 
Nkp46,  F4/80,  CDllb,  Grl,  Ly6G,  CDllc,  and  FoxP3). 

6.  Cancer  stem  cell  content:  flow  cytometric  analyses  with  CD133, 
CD44,  CD24,  and  CD  117. 

7.  Angiogenesis:  CD31  and  CD34. 

8.  Apoptosis,  DNA  damage:  ApopTag,  CC3  positivity,  PARP 
cleavage,  or  histone  H2AX  phosphorylation. 

9.  Toxicity:  histological  analysis  of  liver,  lung,  and  kidney  injury 

10.  TGFP  signaling:  immunodetection  of  phosphorylated  Smad2/3. 

3-30 

10% 

Subtask  5 

Analyze  data  using  statistical  methods;  replicate  experiments  if 
necessary,  prepare  and  submit  manuscripts. 

3-36 

5% 

Milestone  Achieved 

Verified  therapeutic  efficacy  of  CTGF,  CTSK,  FN1,  and  LOXL2 
inhibitors. 

32 

5% 

Specific  Aim  2  (specified  in  proposal) 

Timeline 

Site  1 

Major  Task  1 

Determine  the  effect  of  COL11A1  knockdown  in  human 
carcinoma-associated  fibroblasts 

Months 

Cedars-Sinai 

Medical 

Center 

Subtask  1 

Amend  approved  IACUC  protocol  5318  (Mouse  Models  of  Tumor 
Microenvironment,  PI:  Orsulic)  for  local  approval  and  send  related 
material  for  DoD’s  approval. 

Upon  award 
notice 

100% 

Subtask  2 

Knock  out  COL11A1  in  human  carcinoma-associated  fibroblasts  using 
CRISPR. 

1-3 

10% 

Subtask  3 

Co-culture  COL11A1  knockout  carcinoma-associated  fibroblasts  with 

3-12 

10% 

2 


ovarian  cancer  cells  under  kidney  capsule  of  nude  mice;  measure  cell 
proliferation,  cell  death  and  other  parameters. 

Subtask  4 

Analyze  data  using  statistical  methods;  replicate  experiments  if 
necessary. 

12-18 

10% 

Milestone  Achieved 

Verified  whether  COL11A1  in  carcinoma-associated  fibroblasts  is 
essential  for  the  tumor  promoting  effects  in  a  paracrine  manner. 

18 

10% 

Major  Task  2 

Determine  the  potential  of  COL11A1  as  a  therapeutic  target 

Months 

Cedars-Sinai 

Medical 

Center 

Subtask  1 

Amend  approved  IACUC  protocol  5318  (Mouse  Models  of  Tumor 
Microenvironment,  PI:  Orsulic)  for  local  approval  and  send  related 
material  for  DoD’s  approval. 

Upon  award 
notice 

100% 

Subtask  2 

Purchase  FVB  mice,  drugs  and  reagents;  plan  experiments. 

1-2 

50% 

Subtask  3 

Implant  FVB  mice  with  mouse  ovarian  cancer  cells. 

2-25 

50% 

Subtask  4 

Treat  mice  with  COL11A1  neutralizing  antibody. 

Assess  therapeutic  efficacy  as  in  Aim  1,  Task  4. 

3-30 

50% 

Subtask  5 

Analyze  data  using  statistical  methods;  replicate  experiments  if 
necessary;  prepare  and  submit  manuscripts. 

3-36 

10% 

Milestone  Achieved 

Verified  whether  COL11A1  is  promising  as  a  therapeutic  target  with 
high  specificity  for  activated  carcinoma-associated  fibroblasts. 

36 

10% 

Specific  Aim  3  (specified  in  proposal) 

Timeline 

Site  1 

Major  Task  1 

Assess  the  effect  of  differentiating  carcinoma-associated 
fibroblasts  into  cartilage  on  tumor  progression  and 
chemosensitivity 

Months 

Cedars-Sinai 

Medical 

Center 

Subtask  1 

Amend  approved  IACUC  protocol  5318  (Mouse  Models  of  Tumor 
Microenvironment,  PI:  Orsulic)  for  local  approval  and  send  related 
material  for  DoD’s  approval. 

Upon  award 
notice 

100% 

Subtask  2 

Purchase  FVB  mice,  drugs  and  reagents;  plan  experiments. 

1-2 

30% 

Subtask  3 

Implant  FVB  mice  with  mouse  ovarian  cancer  cells. 

2-25 

30% 

3 


Subtask  4 

Treat  mice  with  recombinant  collagen  II,  rAAV-FLAG-Sox9,  and 
dexamethasone. 

Assess  therapeutic  efficacy  as  in  Aim  1,  Task  4. 

3-30 

30% 

Subtask  5 

Analyze  data  using  statistical  methods;  replicate  experiments  if 
necessary;  prepare  and  submit  manuscripts. 

4-36 

10% 

Milestone  Achieved 

Verified  whether  agents  that  induce  terminal  differentiation  of 
activated  carcinoma-associated  fibroblasts  are  effective  in  attenuating 
tumor  growth  and  increasing  chemosensitivity. 

36 

10% 

What  was  accomplished  under  these  goals? 


1 )  major  activities 

Our  activities  focused  on  obtaining  the  required  regulatory  documents  for  work  with  human 
tissues  and  mouse  models,  planning  experiments,  purchasing  reagents,  and  conducting 
experiments.  Due  to  the  premature  departure  of  the  postdoctoral  fellow  who  conducted  the 
preclinical  tests  in  mice,  we  have  not  yet  analyzed  the  majority  of  the  tissues  that  have  been 
collected  from  these  experiments. 

2)  specific  objectives 

We  have  fully  optimized  the  preclinical  trial  design  in  our  immunocompetent  mouse  model  of 
ovarian  cancer.  We  also  developed  a  rapid  and  efficient  pipeline  for  the  collection  of  ascites  and 
tissues  from  tumors  and  major  organs  that  may  be  affected  by  anti-fibrotic  agents  (liver,  spleen, 
pancreas,  diaphragm,  kidney,  lung,  heart/aorta,  lymph  nodes,  Peyer’s  patches,  brain, 
ovary/oviduct/uterus  and  bone  marrow)  and  processing  for  immunohistochemistry, 
immunofluorescence,  and  RNA  analysis. 

3 )  significant  results  or  key  outcomes,  including  major  findings,  developments,  or  conclusions 
AIM  1.  Improve  therapeutic  efficacy  by  targeting  processes  involved  in  CAF  activation 

1.  Prioritization  of  therapeutic  targets.  We  prioritized  therapeutic  targets  that  are  likely  to  have 
high  specificity  for  activated  CAFs  and  minimal  toxicity  in  normal  tissues.  The  results  of  our  in 
silico  studies  have  been  published  with  acknowledgment  of  the  grant  support  (Jia  et  al.,  Cancer 
Letters,  2016).  This  study  revealed  that  very  few  targets  are  truly  CAF-specific  and  that  thorough 
validation  in  human  cancer  and  normal  tissues  will  be  necessary  to  identify  the  most  effective 
targets.  We  are  currently  validating  expression  of  our  selected  targets  by  immunohistochemistry 
and  in  situ  hybridization  in  human  primary,  metastatic,  and  recurrent  ovarian  cancers  as  well  as 
normal  tissues  (normal  ovary,  fallopian  tube,  lymph  node,  tonsil,  and  kidney). 

2.  Preclinical  trial  design.  In  order  to  optimize  the  preclinical  protocol  of  cisplatin  in 
combination  with  anti-fibrotic  agents  in  our  immunocompetent  mouse  model  of  ovarian  cancer, 
we  studied  the  efficacy  of  cisplatin  in  combination  with  the  Rho-associated,  coiled-coil 
containing  protein  kinase  (ROCK)  inhibitor  Y-27632.  Studies  have  shown  that  Y-27632:  1) 
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enhances  the  efficacy  of  cisplatin-induced  toxicity  in  ovarian  cancer  cell  lines  in  vitro  [1,  2],  yet 
reduces  neurotoxicity  in  a  mouse  model  of  cisplatin-induced  peripheral  neuropathy  [3];  2) 
inhibits  myofibroblast  contraction  in  collagen  contraction  assays  [4-9];  3)  prevents  fibroblast 
proliferation  in  primary  cancer  cell  culture  (Georgetown  method  [10]);  4)  reduces  fibrosis  in 
mouse  models  [11-14];  and  5)  decreases  cancer  growth  and  metastasis  in  various  cancer  models 
[15-24]  although  these  studies  did  not  address  whether  the  effects  of  Y-27632  were  primarily  due 
to  changes  in  cancer  cells  or  CAFs.  During  the  course  of  our  experiments,  three  groups  working 
on  mouse  models  of  pancreatic  cancer  have  published  results  demonstrating  that  ROCK 
inhibition  delays  cancer  progression  due  to  inactivation  of  CAFs  [25-27].  Currently,  several 
ROCK  inhibitors  are  in  clinical  trials  for  vascular  diseases  (clinicaltrials.gov). 

The  goal  of  our  study  was  to  test  the  effects  of  ROCK  inhibition  with  or  without  cisplatin  using 
the  following  parameters: 

1 .  tumor  growth  and  ascites  accumulation  (survival  endpoint) 

2.  alterations  in  the  tumor  microenvironment  (fibrosis,  stem  cell  enrichment,  enrichment  in 
specific  immune  cell  types) 

3.  alterations  in  normal  organs  (liver,  spleen,  pancreas,  diaphragm,  kidney,  lung,  heart/aorta, 
lymph  nodes,  Peyer’s  patches,  brain,  ovary/oviduct/uterus  and  bone  marrow) 

Experiment  1.  Short-term  treatment  to  inhibit  ovarian  cancer  formation  (7  treatments,  26  days; 

20  mice,  4  groups).  The  goal  of  this  experiment  was  to  assess  the  effect  of  short-term  treatment 
on  ovarian  cancer  progression  before  it  reached  the  advanced  stage.  FVB  mice  were  inoculated 
with  100  pi  of  5xl05  syngeneic  Br-Luc  ovarian  cancer  cells  by  intraperitoneal  (i.p.)  injection. 
Twelve  days  after  cancer  cell  injection,  mice  were  randomized  and  treated  i.p.  with  vehicle 
(PBS),  Y-27632  (10  mg/kg),  cisplatin  (1  mg/kg),  or  a  combination  of  Y-27632  (10  mg/kg)  and 
cisplatin  (1  mg/kg),  in  a  total  volume  of  100  pi,  3x/weekly  over  26  days.  The  results  showed  no 
statistical  difference  in  cancer  onset  between  groups  treated  with  vehicle  or  the  ROCK  inhibitor. 
There  was  also  no  statistical  difference  in  cancer  onset  between  groups  treated  with  the  ROCK 
inhibitor  alone  or  with  a  combination  of  the  ROCK  inhibitor  and  cisplatin. 

Experiment  2.  Long-term  treatment  to  inhibit  ovarian  cancer  formation  (18  treatments,  55  days; 
19  mice,  4  groups).  The  goal  of  this  experiment  was  to  assess  the  effect  of  long-term  treatment 
on  ovarian  cancer  progression  before  it  reached  the 
advanced  stage.  FVB  mice  were  inoculated  with  100 
pi  of  5x10s  syngeneic  Br-Luc  ovarian  cancer  cells  by 
i.p.  injection.  Twelve  days  after  cancer  cell  injection, 
mice  were  randomized  and  treated  i.p.  with  vehicle 
(PBS),  Y-27632  (10  mg/kg),  cisplatin  (1  mg/kg),  or  a 
combination  of  Y-27632  (10  mg/kg)  and  cisplatin  (1 
mg/kg),  in  a  total  volume  of  100  pi,  3x/weekly  over 
55  days.  The  treatment  schedule  of  the  combination 
group  was  reduced  to  lx/weekly  after  4  weeks  of 
treatment  due  to  possible  toxicity  (mice  displayed 
jumpiness).  Kaplan-Meier  survival  curve  is  shown  in 
Fig.  1. 
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Experiment  3.  Treatment  to  reduce  advanced-stage  ovarian  cancer  (2-7  treatments;  2-7  days;  20 
mice,  4  groups).  The  goal  of  this  experiment  was  to  assess  the  effect  of  treatment  on  advanced- 
stage  ovarian  cancer.  FVB  mice  were  inoculated  with  100  pi  of  5xl05  syngeneic  Br-Luc  ovarian 
cancer  cells  by  i.p.  injection.  Mice  with  ascites  (symptom  of  advanced- stage  ovarian  cancer) 
were  treated  i.p.  with  vehicle  (PBS),  Y-27632  (10  mg/kg),  cisplatin  (1  mg/kg),  or  a  combination 
of  Y-27632  (10  mg/kg)  and  cisplatin  ( 1  mg/kg),  in  a  total  volume  of  100  pi.  The  experiment  was 
terminated  due  to  the  early  death  of  mice  that  did  not  receive  cisplatin  in  their  treatment. 

Experiment  4.  Pre-treatment  with  the  ROCK  inhibitor  alone  during  cancer  development  and 
subsequent  treatment  with  the  ROCK  inhibitor  in  combination  with  a  low  dose  of  cisplatin  upon 
advanced- stage  ovarian  cancer  development,  (multiple  treatments,  10  mice,  2  groups).  The  goals 
of  this  experiment  were  to  assess  the  effect  of  long-term  ROCK  inhibition  on  the  onset  of 
advanced-stage  cancer  (accumulation  of  ascites)  and  to  assess  the  effect  of  long-term 
pretreatment  with  the  ROCK  inhibitor  on  tumor  burden  reduction  when  advanced-stage  ovarian 
cancer  is  treated  with  a  low  dose  of  cisplatin.  FVB  mice  were  inoculated  with  100  pi  of  5xl05 
syngeneic  Br-Luc  ovarian  cancer  cells  by  i.p.  injection.  Two  days  after  cancer  cell  inoculation, 
mice  were  treated  with  i.p.  with  vehicle  (PBS)  or  Y-27632  (10  mg/kg)  until  they  developed 
ascites  (symptom  of  advanced-stage  ovarian  cancer).  At  that  point,  mice  were  treated  i.p.  with 
vehicle  (PBS)  +  cisplatin  (0.2  mg/kg)  or  Y-27632  (10  mg/kg)  +  cisplatin  (0.2  mg/kg),  in  a  total 
volume  of  100  pi.  One  mouse  in  the  Y-27632  pretreatment  group  did  not  develop  cancer  by  the 
end  of  the  experiment.  The  results  are  shown  in  Fig.  2. 
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Aim  2.  Increase  specificity  of  targeting  activated  CAFs 

Our  analyses  identified  COL1 1 A1  as  the  most  specific  target  for  activated  CAFs.  We 
hypothesized  that  targeting  COL11A1  function  will  disable  activated  CAFs  with  a  minimal 
effect  on  normal  fibroblasts.  We  used  two  approaches  to  inactivate  COL11A1. 

1.  Knockout  of  COL11A1  in  CAFs.  To  test  whether  78  IT  CAFs  support  growth  of  ovarian 
cancer  cells,  human  ovarian  cancer  cell  line  SKOV3-GFP  was  implanted  subcutaneously  alone 
or  with  781T  CAFs  (10  mice,  5  mice  per  group,  2  flanks  per  mouse).  SKOV3-GFP  cells  co¬ 
injected  with  78 IT  CAFs  formed  larger  tumors  (Fig.  3A)  that  also  grew  well  when  explanted  in 
culture  (Fig.  3B).  Immunohistochemical  analysis  of  the  subcutaneous  tumors  showed  abundant 
lymphocytic  infiltrates  associated  with  SKOV3-GFP  cells  co-injected  with  78 IT  CAFs  but  not 
with  SKOV3-GFP  cells  injected  alone  (data  not  shown),  indicating  that  CAFs  are  not  only 
effective  in  supporting  cancer  cell  growth  but  also  play  a  role  in  the  recruitment  of  other 


|  Time  of  onset  of  ascites 
|  Time  to  morbidty  after  treatment 


Fig.  2.  A  bar  graph  showing 
time  to  disease  onset  during 
pretreatment  with  PBS  or  Y- 
27632  (blue)  and  time  to 
morbidity  during  treatment  with 
cisplatin  or  a  combination  of 
cisplatin  and  Y-27632  (orange). 
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components  of  the  stromal  microenvironment.  Despite  multiple  efforts,  we  were  unable  to 
effectively  knock  down  COL11A1  in  78 IT  CAFs.  We  will  continue  to  test  different  approaches 
to  knocking  down  COL1 1A1  in  78  IT  CAFs  or  in  a  different  line  of  CAFs. 


SKOV3-GFP 
+  781TCAF 


9x7  mm  9x5  mm 
133  mg  124  mg  514  mg  518  mg 


9  # 

13  x  12  mm  12  x  10  mm 


SKOV3-GFP 


SKOV3-GFP 


SKOV3-GFP 
+  781TCAF 


Fig.  3.  CAFs  support  subcutaneous  growth  of  SKIO- 
GFP  human  ovarian  cancer  cells.  A)  Subcutaneous 
tumors  isolated  from  mice  40  days  after  injection  of 
10.000  SKOV3-GFP  cells  with  or  without  20,000 
78 IT  CAFs.  B)  Tumor  explants  in  culture  imaged 
under  GFP  or  visible  light.  Explants  from  SKOV3- 
GFP  tumors  grow  poorly  in  culture,  while  explants 
from  SKOV3-GFP  +  78 IT  CAF  tumors  grow  well  in 
culture  and  consist  of  both  GFP-positive  epithelial 
cells  and  GFP-negative  fibroblasts  (circled). 


2.  Neutralizing  endogenous  COL11A1  with  a  COL1  lAl-specific  antibody.  In  an  experiment 
with  10  mice  treated  with  control  IgG  or  COL11A1  antibody  (2.3mg/ml,  daily  for  10  days),  we 
did  not  see  a  statistically  significant  difference  in  tumor  growth  between  the  two  groups  of  mice, 
indicating  that  this  antibody  (or  this  concentration)  is  not  effective  in  inhibiting  tumor  growth 
(data  not  shown). 

Aim  3.  Induce  activated  CAF-to-cartilage  differentiation. 

We  are  attempting  to  drive  activated  CAFs  into  an  evolutionary  dead-end  by  inducing  terminal 
cell  differentiation.  Bioengineers  have  observed  that  cartilage  cells  spontaneously  de¬ 
differentiate  into  myofibroblasts  upon  growth  on  plastic,  indicating  that  reverse  differentiation 
should  be  possible.  We  identified  that  CAFs  are  poised  for  cartilage  differentiation  if  inhibitory 
differentiation  signals  are  disabled.  To  date,  we  have  tested  dexamethasone  as  a  differentiation 
agent.  We  did  not  see  a  significant  difference  in  tumor  growth  (data  not  shown).  However,  the 
collected  tissues  have  yet  to  be  analyzed  for  cartilage/bone  formation  markers  to  assess  the  level 
of  differentiation. 
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4)  other  achievements 
Nothing  to  report. 

a.  What  opportunities  for  training  and  professional  development  has  the  project 
provided? 
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Nothing  to  report. 

b.  How  were  the  results  disseminated  to  communities  of  interest? 

Nothing  to  report. 

c.  What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

For  Aims  1  and  3,  tumor  and  normal  tissues  (liver,  spleen,  pancreas,  diaphragm,  kidney,  lung, 
heart/aorta,  lymph  nodes,  Peyer’s  patches,  brain,  ovary/oviduct/uterus  and  bone  marrow)  have 
been  collected  from  all  experiments  for  immunohistochemical  (formalin-fixed  and  OCT-frozen) 
and  RNA  (snap-frozen)  analyses  of  markers  of  toxicity,  fibrosis,  stem  cell  enrichment,  and 
enrichment  in  specific  immune  cell  types.  To  understand  the  possible  role  of  Y-27632  in 
increasing  the  accessibility  of  cisplatin,  we  will  conduct  thorough  cell  apoptosis  evaluation  in 
tumors  and  normal  organs.  These  analyses  are  currently  in  process.  The  results  of  Y-27632 
experiments  will  be  validated  and  published.  One  or  two  additional  anti-fibrotic  agents  (Aim  1) 
and  differentiation  agents  (Aim  3)  will  be  tested  during  the  next  reporting  period. 

For  Aim  2,  we  will  try  different  approaches  to  knocking  down  COL11A1  in  CAFs  to  test  the 
function  of  COL11A1  in  the  promotion  of  tumor  growth  by  CAFs.  We  will  also  try  different 
concentrations  and  routes  of  delivery  for  the  COL11A1  antibody  (10X  concentration  and  direct 
intratumoral  delivery  into  a  localized  subcutaneous  tumor)  to  rule  out  the  possibility  that 
COL11A1  did  not  show  an  effect  on  tumor  growth  because  of  a  low  concentration  and/or  rapid 
degradation  upon  intraperitoneal  delivery. 

3.  IMPACT: 

a.  What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the 
project? 

Nothing  to  report. 

b.  What  was  the  impact  on  other  disciplines? 

Nothing  to  report. 

c.  What  was  the  impact  on  technology  transfer? 

Nothing  to  report. 

d.  What  was  the  impact  on  society  beyond  science  and  technology? 

Nothing  to  report. 
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4.  CHANGES/PROBLEMS: 


a.  Changes  in  approach  and  reasons  for  change.  No 

b.  Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them. 

We  successfully  resolved  an  unanticipated  challenge  during  the  project  period.  We  were  pleased 
to  learn  that  Dr.  Dongyu  Jia,  a  funded  postdoctoral  fellow  on  the  project,  was  named  to  a  faculty 
position  in  the  University  System  of  Georgia.  We  were,  with  careful  planning  and  support  from 
our  Human  Resources  department,  able  to  post  and  hire  a  highly  skilled  postdoctoral  fellow  with 
<  2  month  delay.  Dr.  Marcela  Haro  worked  with  Dr.  Orsulic  and  the  Research  Associate 
assigned  to  the  project  to  assume  the  duties  of  the  study.  There  was  minimal  disruption  to  the 
project. 

c.  Changes  that  had  a  significant  impact  on  expenditures. 

The  changes  detailed  in  section  above  had  no  significant  impact  on  overall  expenditures. 

d.  Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals, 
biohazards,  and/or  select  agents.  No 

e.  Significant  changes  in  use  or  care  of  human  subjects.  No 

f.  Significant  changes  in  use  or  care  of  vertebrate  animals.  No 

g.  Significant  changes  in  use  of  biohazards  and/or  select  agents.  No 

5.  PRODUCTS: 

Nothing  to  report. 

a.  Publications,  conference  papers,  and  presentations 

i.  Journal  publications. 

Jia  D,  Liu  Z,  Deng  N,  Tan  TZ,  Huang  R  Y-J,  Taylor-Harding  B,  Cheon  DJ,  Lawrenson  K, 
Wiedemeyer  WR,  Walts  AE,  Kalian  BY,  Orsulic  S.  A  COL1 1  Al-correlated  pan-cancer  gene 
signature  of  activated  fibroblasts  for  the  prioritization  of  therapeutic  targets.  Cancer  Letters 
2016;  382:203-214.  (Published,  acknowledged  grant  funding) 

ii.  Books  or  other  non-periodical,  one-time  publications.  N/A 

iii.  Other  publications,  conference  papers,  and  presentations. 

Invited  talks  (acknowledged  grant  funding) 

Sandra  Orsulic:  Tumor  Microenvironment.  Molecular  and  Cellular  Biology  Wednesday  Seminar 
Series,  Baylor  College  of  Medicine.  Houston,  TX.  January  6,  2016. 

Sandra  Orsulic:  Signatures  of  Stromal  Activation  in  Cancer.  Molecular  Pathology  Seminar 
Series.  Johns  Hopkins  University  School  of  Medicine.  Baltimore,  MD.  November  16,  2016. 
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b.  Website(s)  or  other  Internet  site(s).  N/A 

c.  Technologies  or  techniques.  N/A 

d.  Inventions,  patent  applications,  and/or  licenses.  N/A 

e.  Other  Products.  N/A 

6.  PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 

a.  What  individuals  have  worked  on  the  project? 


Name: 

Sandra  Orsulic 

Project  Role: 

PI 

Nearest  person  month  worked: 

2.00 

Contribution  to  Project: 

Dr.  Orsulic  oversaw  projects  for  all  three  specific  aims,  including 
experimental  design,  execution,  and  data  analysis  and  interpretation.  She 
wrote  the  manuscript  (Jia  et  al,  Cancer  Letters,  2017)  and  prepared 
presentations  as  well  as  the  progress  report. 

Name: 

Beth  Karlan 

Project  Role: 

Collaborator 

Nearest  person  month  worked: 

0.12 

Contribution  to  Project: 

Dr.  Karlan  advised  on  the  translational  aspects  of  the  proposal  and 
participated  in  experimental  design. 

Name: 

Dongyu  Jia,  PhD 

Project  Role: 

Postdoctoral  Fellow 

Nearest  person  month  worked: 

0.60 

Contribution  to  Project: 

Dr.  Jia  conducted  all  experiments  that  involved  testing  different 
combinations  of  treatments  in  the  immunocompetent  mouse  model  of 
ovarian  cancer  and  assisted  in  data  acquisition,  analysis,  and 
interpretation  as  well  as  in  the  writing  of  the  published  manuscript. 

Name: 

Marcela  Haro,  PhD 

Project  Role: 

Postdoctoral  Fellow 

Nearest  person  month  worked: 

4.50 

Contribution  to  Project: 

Dr.  Haro  conducted  all  experiments  that  involved  testing  different 
combinations  of  treatments  in  the  immunocompetent  mouse  model  of 
ovarian  cancer  and  assisted  in  data  acquisition  and  analysis. 
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Name: 

Barbie  Taylor-Harding,  PhD 

Project  Role: 

Research  Associate 

Nearest  person  month  worked: 

3.60 

Contribution  to  Project: 

Dr.  Taylor-Harding  conducted  all  experiments  that  involved  testing 
different  combinations  of  treatments  in  the  immunocompetent  mouse 
model  of  ovarian  cancer  and  assisted  in  data  acquisition  and  analysis. 

b.  Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or  senior/key 
personnel  since  the  last  reporting  period?  No 

c.  What  other  organizations  were  involved  as  partners?  None 

7.  SPECIAL  REPORTING  REQUIREMENTS 

a.  COLLABORATIVE  AWARDS:  N/A 

b.  QUAD  CHARTS:  N/A 

8.  APPENDICES: 

Reprint  of  the  manuscript  Jia  et  al,  Cancer  Letters,  2016 
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Although  cancer-associated  fibroblasts  (CAFs)  are  viewed  as  a  promising  therapeutic  target,  the  design 
of  rational  therapy  has  been  hampered  by  two  key  obstacles.  First,  attempts  to  ablate  CAFs  have  re¬ 
sulted  in  significant  toxicity  because  currently  used  biomarkers  cannot  effectively  distinguish  activated 
CAFs  from  non-cancer  associated  fibroblasts  and  mesenchymal  progenitor  cells.  Second,  it  is  unclear  whether 
CAFs  in  different  organs  have  different  molecular  and  functional  properties  that  necessitate  organ- 
specific  therapeutic  designs.  Our  analyses  uncovered  COL1 1 A1  as  a  highly  specific  biomarker  of  activated 
CAFs.  Using  COL11A1  as  a  ‘seed’,  we  identified  co-expressed  genes  in  13  types  of  primary  carcinoma  in 
The  Cancer  Genome  Atlas.  We  demonstrated  that  a  molecular  signature  of  activated  CAFs  is  conserved 
in  epithelial  cancers  regardless  of  organ  site  and  transforming  events  within  cancer  cells,  suggesting  that 
targeting  fibroblast  activation  should  be  effective  in  multiple  cancers.  We  prioritized  several  potential 
pan-cancer  therapeutic  targets  that  are  likely  to  have  high  specificity  for  activated  CAFs  and  minimal  tox¬ 
icity  in  normal  tissues. 

©  2016  The  Author(s).  Published  by  Elsevier  Ireland  Ltd.  This  is  an  open  access  article  under  the  CC  BY- 

NC-ND  license  (http://creativecommons.Org/licenses/by-nc-nd/4.0/). 


Introduction 

Under  normal  physiological  conditions,  collagen-rich  fibro¬ 
blasts  maintain  tissue  architecture  and  serve  as  a  barrier  to  epithelial 
cell  migration.  However,  cancer  cells  have  the  ability  to  convert 
the  surrounding  fibroblasts  into  activated  CAFs,  which  secrete  spe¬ 
cific  collagens,  growth  factors,  and  enzymes  that  promote  cancer 
growth,  angiogenesis,  invasion,  and  metastasis  [1-3].  At  the  same 
time,  these  CAFs  suppress  anticancer  immunity,  confer  drug  resis¬ 
tance  and/or  limit  the  access  of  chemotherapies,  anti-angiogenic 
therapies,  and  immunotherapies  [  1  — 3 ].  Although  the  exact  mecha¬ 


Abbreviations:  aSMA,  a-smooth  muscle  actin;  CAFs,  cancer-associated  fibro¬ 
blasts;  CSPG4,  chondroitin  sulfate  proteoglycan  4;  ECM,  extracellular  matrix;  LOX, 
lysyl  oxidase;  EMT,  epithelial-mesenchymal  transition;  FAP,  fibroblast  activation 
protein;  PALLD,  palladin;  PDGFRa,  platelet-derived  growth  factor  receptor  a;  PDPN, 
podoplanin;  TGFp,  transforming  growth  factor  (1;  TNC,  tenascin-C;  PRECOG,  PREdiction 
of  Clinical  Outcomes  from  Genomic  Profiles. 
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nisms  by  which  activated  CAFs  contribute  to  such  diverse  aspects 
of  cancer  progression  are  unclear,  it  is  thought  that  fibroblasts 
together  with  increased  collagen  deposition  and  altered  extracel¬ 
lular  matrix  (ECM)  remodeling  serve  as  a  rich  depot  of  cancer- 
promoting  growth  factors,  cytokines,  and  chemokines  [1,2], 
Additionally,  altered  levels  of  enzymes  responsible  for  collagen 
cross-link  formation,  such  as  lysyl  oxidase  (LOX)  [4],  increase  tissue 
stiffness  and  modify  mechanotransduction  resulting  in  the  reorga¬ 
nization  of  loose  connective  tissue  into  tense  linear  tracks  of  fibers 
that  serve  as  highways  to  promote  chemotaxis  of  cancer  cells 
[5,6], 

Recognizing  the  crucial  role  of  CAFs  in  most  aspects  of  cancer 
progression,  it  has  been  proposed  that  rational  anticancer  therapy 
design  should  not  only  target  the  cancer  cells  but  also  the  CAFs  [7,8  j. 
Unlike  cancer  cells,  CAFs  are  genetically  stable  [9],  which  reduce  the 
risk  of  therapy-induced  clonal  selection,  resistance,  and  cancer  re¬ 
currence.  Furthermore,  targeting  CAFs  could  potentially  affect 
multiple  biochemical  pathways  to  prevent  cancer  progression 
and  recurrence.  CAF-targeting  therapeutic  approaches  in  different 
experimental  mouse  cancer  models  have  been  shown  to  improve 


http://dx.doi. org/10.1 01 6/j.canlet.201 6.09.001 

0304-3835/©  2016  The  Author(s).  Published  by  Elsevier  Ireland  Ltd.  This  is  an  open  access  article  under  the  CC  BY-NC-ND  license  (http://creativecommons.org/licenses/ 
by-nc-nd/4.0/). 


204 


D.  Jia  et  al. /Cancer  Letters  382  (2016)  203-214 


tumoral  immune  response,  intratumoral  drug  delivery,  and  thera¬ 
peutic  efficacy  [10-15],  These  studies  confirm  the  key  role  of  CAFs 
in  cancer  progression  and  demonstrate  their  effectiveness  as  a  ther¬ 
apeutic  target.  However,  targeting  CAFs  in  some  cancer  models 
actually  promoted  cancer  progression.  For  example,  depletion  of 
aSMA+  stroma  in  a  mouse  pancreatic  cancer  model  resulted  in  in¬ 
creased  cancer  aggressiveness,  enhanced  hypoxia  and  epithelial- 
mesenchymal  transition  (EMT),  suppressed  anticancer  immunity, 
and  reduced  survival  [16], 

The  contradictory  results  in  different  cancer  models  could  be  ex¬ 
plained  by  different  roles  of  CAFs  in  different  cancer  types,  i.e.  CAFs 
could  be  promoting  breast  cancer  and  inhibiting  pancreatic  cancer. 
Alternatively,  in  all  cancer  types  CAFs  prevent  cancer  progression 
until  they  receive  activating  signals  from  cancer  cells  and  convert 
into  ‘activated  CAFs’,  which  in  turn  confer  invasive  and  metastatic 
abilities  upon  cancer  cells  [7].  Therapies  that  target  all  CAFs  are  coun¬ 
terproductive  and  likely  to  result  in  the  death  of  normal  fibroblasts 
and  significant  toxicity.  Preferential  targeting  of  activated  CAFs  has 
been  challenging  because  activated  CAFs  are  poorly  understood  at 
the  molecular  level.  During  activation,  CAFs  exhibit  phenotypic 
changes  that  partially  overlap  with  myofibroblastic  changes  during 
wound  healing,  inflammation,  and  fibrosis,  including  secretion  of 
specific  ECM  components,  cytokines  and  growth  factors  [1,17], 
Several  markers  have  been  used  to  distinguish  activated  from  non- 
activated  CAFs:  a-smooth  muscle  actin  (aSMA,  encoded  by  gene 
ACTA2)  [3],  fibroblast  activation  protein  (FAP)  [12],  podoplanin 
(PDPN)  [18],  palladin  (PALLD)  [19,20],  tenascin-C  (TNC)  [21  ],  platelet- 
derived  growth  factor  receptor  a  (PDGFRa)  [22,23],  and  chondroitin 
sulfate  proteoglycan  4  (CSPG4)  [24].  However,  these  markers  are  fre¬ 
quently  expressed  in  other  cells  within  the  cancer  stroma,  such  as 
vascular  smooth  muscle  cells,  pericytes,  and  mesenchymal  stem  cells. 
This  lack  of  specificity  could  pose  problems  in  therapeutic  target¬ 
ing  and  underscores  the  need  to  better  understand  the  molecular 
characteristics  of  activated  CAFs  in  order  to  develop  more  precise 
and  less  toxic  targeted  therapies. 

COL11A1  encodes  the  al  chain  of  collagen  XI,  a  minor  fibrillar 
collagen  expressed  by  chondrocytes  and  osteoblasts  but  not  qui¬ 
escent  fibroblasts  [25,26],  The  absence  of  functional  collagen  XI  leads 
to  abnormally  thickened  cartilage  and  tendon  fibrils,  suggesting  the 
role  of  collagen  XI  in  maintaining  proper  fibril  diameter  [27,28], 
Studies  have  demonstrated  that  COLIIAI  mRNA  is  markedly  el¬ 
evated  in  cancers  of  the  oral  cavity/pharynx,  head  and  neck,  breast, 
lung,  esophagus,  stomach,  pancreas,  colon,  and  ovary,  but  not  in 
matched  normal  tissues  (reviewed  in  [25,26]).  COLIIAI  has  been 
identified  as  part  of  gene  signatures  associated  with  adverse  clin¬ 
ical  outcomes  including  resistance  to  neoadjuvant  therapy  in  breast 
cancer  [29],  time  to  recurrence  in  glioblastoma  [30],  poor  survival 
in  kidney  cancer  [31],  and  time  to  recurrence  and  overall  survival 
in  ovarian  cancer  [32,33].  In  situ  hybridization  in  ovarian  cancer 
and  immunohistochemistry  in  pancreatic  cancer  revealed  that 
COL11A1  mRNA  and  pro-protein  are  primarily  expressed  in  CAFs 
[25,32],  The  restricted  expression  of  COLIIAI  in  normal  tissues  and 
its  enrichment  in  CAFs  during  cancer  progression  combined  with 
its  association  with  adverse  clinical  outcomes  in  multiple  types  of 
cancer  support  its  candidacy  as  a  specific  marker  of  fibroblast  ac¬ 
tivation  in  diverse  cancers.  Here,  we  explore  the  suitability  of  COLIIAI 
as  a  pan-cancer  marker  of  activated  CAFs  and  use  it  as  a  ‘seed’  to 
identify  the  transcription  signature  of  activated  CAFs  in  13  epithe¬ 
lial  cancer  types.  We  show  that  the  COLIIAI -coexpressed  gene  set 
is  highly  conserved  in  these  13  cancer  types,  indicating  that  the  fi¬ 
broblast  reaction  to  cancer  cells  is  independent  of  the  organ  site- 
of-origin  and  of  the  transforming  events  within  cancer  cells.  Finally, 
by  combining  drug  target  databases  with  cancer  vs.  normal  tissue 
expression  databases,  we  identify  several  potential  therapeutic 
targets  that  should  have  high  specificity  for  activated  CAFs  and 
minimal  toxicity  in  normal  tissues. 


Materials  and  methods 

Human  tissues 

Studies  involving  human  tissue  samples  were  approved  by  the  Cedars-Sinai  In¬ 
stitutional  Review  Board  (IRB  15425).  The  samples  included  a  tissue  microarray  from 
42  patients  with  matched  primary,  metastatic,  and  recurrent  ovarian  cancer. 

In  situ  hybridization 

The  RNA  hybridization  kit  (RNAscope  2.0  FFPE  Assay)  and  probes  for  COLIIAI, 
the  bacterial  gene  dapB  (negative  control),  and  the  housekeeping  gene  HPRT  (pos¬ 
itive  control),  were  from  Advanced  Cell  Diagnostics,  Inc.  Formalin-fixed,  paraffin- 
embedded  tissue  section  slides  were  processed  by  the  Cedars-Sinai  Biobank  and 
Translational  Research  Core  following  the  protocol  provided  with  the  RNAscope  In 
Situ  Hybridization  kit  from  Advanced  Cell  Diagnostics,  Inc.  The  slides  were  coun- 
terstained  with  Mayer’s  hematoxylin. 

Immunohistochemistry 

lmmunohistochemical  detection  of  aSMA  was  performed  on  formalin-fixed, 
paraffin-embedded  tissue  sections  using  the  protocol  provided  with  the  pre¬ 
diluted  asm-1  clone  antibody  from  Leica  Microsystems.  Staining  was  done  by  the 
Cedars-Sinai  Pathology  Service  on  the  Ventana  Benchmark  Ultra  automated  slide 
Stainer.  The  staining  was  visualized  using  the  Ventana  OptiView  DAB  Detection  System. 
The  slides  were  counterstained  with  Mayer's  hematoxylin. 

In  vitro  co-culture  experiments 

The  ovarian  cancer  cell  lines  OVCAR3-GFP,  KURAMOCHI-GFP  and  OVSAHO- 
GFP  were  maintained  in  RPMI-1640  (Corning)  supplemented  with  10%  fetal  bovine 
serum  (FBS)  and  lx  penicillin/streptomycin  (Corning).  Cell  line  authenticity  was  con¬ 
firmed  by  Laragen  using  the  short  tandem  repeat  (STR)  method.  The  immortalized 
normal  ovarian  fibroblasts  INOF-tdTomato  cell  lines  [34]  were  maintained  in  a  1:1 
ratio  of  MCDB  105  (Sigma-Aldrich)  and  Medium  199  (GIBCO)  with  10%  FBS,  50  U/ 
ml  penicillin  and  50  ug/ml  streptomycin.  Immortalized  normal  ovarian  fibroblasts 
and  ovarian  cancer  cells  were  co-cultured  in  1%  FBS  supplemented  media  (1:1:2  ratio 
of  MCDB  105,  Medium  199  and  RPMI-1640)  using  6-well  plates,  either  by  directly 
plating  ovarian  cancer  cells  ( 1 05  cells/well)  onto  a  70%  confluent  layer  of  normal 
ovarian  fibroblasts  or  onto  a  0.4  pm  Transwell  membrane.  Media  were  replaced  every 
2  days.  After  4  days  of  co-culture,  GFP-labeled  ovarian  cancer  cells  and  tdTomato- 
labeled  fibroblasts  were  separated  by  FACS  in  PBS  with  0.5%  FBS.  RNA  extraction  from 
fibroblasts  and  ovarian  cancer  cell  lines  was  performed  using  the  RNeasy  Mini  Kit 
(Qiagen)  and  reverse  transcribed  to  cDNA  using  the  Quantitect  Reverse  Transcrip¬ 
tion  Kit  (Qiagen).  For  qRT-PCR,  50  ng  of  cDNA  was  mixed  with  COLIIAI  primers 
(Forward:  5'-GACTATCCCCTCTTCAGAACTGTTAAC-3';  Reverse:  5'-  CTTCTATCAAGTGG 
1 1TCGTGGT 1 T-3')  and  the  iQSYBR-Green  Supermix  (BioRad)  and  run  on  the  CFX96 
Real-Time  System  (BioRad).  Data  were  analyzed  using  the  2_ACT  method  and  nor¬ 
malized  to  INOF-tdTomato  control  to  present  the  fold  change  ratios.  All  mRNA  data 
were  normalized  to  RPL32  expression  (Forward:  5'-ACAAAGCACATGCTGCCCAGTG- 
3';  Reverse:  5'-TTCCACGATGGCTTTGCGGTTC-3').  The  statistical  analyses  were 
performed  using  GraphPad  Prism  (version  6.0;  GraphPad  Software).  The  unpaired 
t  test  was  used  for  data  analyses. 

Public  database  portals  and  dataset  analyses 

Data  from  public  portals  were  used  as  provided  by  individual  portals  without 
additional  processing  or  normalization,  unless  otherwise  indicated.  Box  plots  of 
COLIIAI  expression  in  normal  tissues  and  cancers  were  generated  using  the  Gene 
Expression  across  Normal  and  Tumor  tissue  (GENT)  portal  (medical- 
genome.kribb.re.kr/GENT)  in  which  data  from  multiple  datasets  were  processed  and 
normalized  as  previously  described  [35],  COLl  1A1  expression  level  diagrams  for  in¬ 
flammatory  bowel  disease,  lung  fibrosis,  and  cancers  of  the  colon  and  lung  were 
generated  using  the  R2  MegaSampler  public  portal  (hgserverl.amc.nl/cgi-bin/r2/ 
main.cgi).  A  description  of  the  methods  used  for  data  processing  and  normalization 
is  available  through  the  portal.  Survival  z-scores  for  individual  genes  and  cancer  types 
were  obtained  from  the  PREdiction  of  Clinical  Outcomes  from  Genomic  Profiles 
(PRECOG)  portal  (precog.stanford.edu).  Methods  for  calculating  PRECOG  z-scores  have 
been  published  [36],  Ranking  of  the  COEllAl-correlated  genes  in  13  TCGA  carcino¬ 
ma  types  was  determined  using  data  from  individual  cancer  datasets  that  were 
processed  by  cBio  Portal  (cbioportal.org)  as  previously  described  [37],  Kaplan- 
Meier  survival  plots  and  plots  of  COLIIAI  expression  in  individual  molecular  subtypes 
of  ovarian  carcinoma  were  generated  using  the  ovarian  cancer  microarray  gene  ex¬ 
pression  database  CSIOVDB  (csibio.nus.edu.sg/CSIOVDB/CSIOVDB.html),  which  has 
been  previously  described  [38],  The  dataset  for  fibroblast  and  ovarian  epithelial  cell 
co-culture  was  imported  from  the  Gene  Expression  Omnibus  (ncbi.nlm.nih.gov/ 
geo).  The  Euclidean  distance  clustering  analysis  heatmap  for  the  e-mtab-991  [39] 
and  GSE40595  [40]  datasets  was  generated  using  the  public  R2  GeneSet  Clustering 
Analysis  portal  (hgserverl.amc.nl/cgi-bin/r2/main.cgi),  which  also  describes  methods 
that  were  used  to  process  and  normalize  data  from  datasets  included  in  the  portal. 
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Results 

COL1 1  Al  is  expressed  in  a  subset  of  aSMA-positive  CAFs  and  can  be 
induced  in  normal  fibroblasts  by  the  presence  of  cancer  cells 

To  determine  if  COLIIAI  expression  is  associated  with  fibro¬ 
blast  activation,  we  used  aSMA  as  a  marker  of  activated  CAFs  [3], 
Comparison  of  aSMA  immunohistochemistry  and  COLIIAI  in  situ 
hybridization  in  a  tissue  microarray  consisting  of  primary,  meta¬ 
static  and  recurrent  ovarian  cancers  from  42  patients  showed  that 
COLIIAI  is  expressed  in  a  subset  of  aSMA+  CAFs  (Fig.  1A).  Unlike 
aSMA,  COLIIAI  was  not  expressed  in  blood  vessels  (red  arrows)  or 
in  fibroblasts  surrounding  the  cancer  (blue  arrows)  (Fig.  1A). 


In  sections  of  metastatic  ovarian  cancer,  we  observed  that 
COL  It  A  3  -positive  cells  are  confined  to  the  intratumoral  and  imme¬ 
diate  peritumoral  CAFs  (Fig.  IB),  suggesting  that  COLiiAi  expression 
may  be  induced  by  cues  received  from  epithelial  cancer  cells.  To  test 
if  cancer  cells  can  induce  COLilAl  expression  in  fibroblasts,  we 
co-cultured  immortalized  normal  ovarian  fibroblasts  (INOFs)  with 
three  different  ovarian  cancer  cell  lines  (OVSAHO,  OVCAR3,  and 
KURAMOCHI).  COLIIAI  expression  in  INOFs  was  most  strongly 
induced  by  direct  co-culture  with  ovarian  cancer  cell  lines  al¬ 
though  weak  induction  occurred  by  indirect  co-culture  on  a 
Transwell  membrane  (Fig.  1C).  The  induction  of  COLIIAI  in  fibro¬ 
blasts  in  the  presence  of  cancer  cells  was  confirmed  by  analysis 
of  the  public  expression  dataset  GSE52104  in  which  two  types  of 
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Fig.  1.  COLIIAI  is  expressed  in  CAFs.  (A)  Comparison  of  aSMA  immunohistochemistry  and  COLIIAI  in  situ  hybridization  in  a  metastatic  ovarian  cancer  sample.  Red  arrows 
indicate  blood  vessels.  Blue  arrows  indicate  fibroblasts  surrounding  the  tumor  nodule.  A  high  magnification  of  peritumoral  and  intratumoral  regions  in  the  black  rectangles 
is  shown  in  panels  on  the  left.  IHC,  immunohistochemistry;  ISH,  in  situ  hybridization;  BV,  blood  vessel;  NC,  necrosis.  (B)  Distribution  of  in  situ  hybridization  COL11A1- 
positive  CAFs  in  relation  to  cancer  cells.  The  estimated  percent  of  COLIIAI  -positive  CAFs  is  shown  on  the  right.  The  image  is  representative  of  metastatic  and  recurrent 
ovarian  cancer  samples,  which  typically  express  higher  levels  of  COLIIAI  than  primary  ovarian  cancers.  (C)  Quantitative  RT-PCR  levels  of  COLIIAI  in  sorted  (FACS)  immor¬ 
talized  normal  ovarian  fibroblasts  (INOFs)  grown  alone  or  co-cultured  with  ovarian  cancer  cell  lines  (OVSAHO,  OVCAR3,  KURAMOCHI)  that  were  either  separated  from  INOFs 
by  a  Transwell  membrane  or  directly  mixed  with  INOFs.  Statistical  analyses  were  performed  between  INOFs  grown  alone  and  INOFs  co-cultured  with  ovarian  cancer  cells 
(*p<0.05;  **p<0.01;  ***p<  0.001;  ns,  not  significant).  Error  bars  indicate  standard  deviation.  (D)  Levels  of  COLIIAI  in  the  GSE52104  expression  dataset  in  which  mesen¬ 
chymal  stem  cells  (MSCs)  or  immortalized  normal  ovarian  fibroblasts  (INOFs)  were  either  cultured  alone  or  co-cultured  with  IOSE4  normal  epithelial  cells  (IOSE)  or  HEYA8 
epithelial  ovarian  cancer  cells  (EOC)  using  a  Transwell  membrane.  Inverse-log2  values  of  the  Robust  Multi-array  Average  (RMA)  scores  from  different  COLIIAI  probes  were 
averaged,  then  log2-transformed.  The  data  were  extracted  for  statistical  analyses  using  GraphPad  Prism  6.  Data  are  represented  as  the  mean  ±  SEM.  Intergroup  differences 
were  assessed  by  the  Student’s  t-test.  *p  <  0.05;  ****p  <  0.0001.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version 
of  this  article.) 
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presumptive  cancer-associated  fibroblast  precursor  cells,  mesen¬ 
chymal  stem  cells  (MSCs)  and  immortalized  normal  ovarian 
fibroblasts  (INOFs),  were  either  cultured  alone  or  co-cultured  with 
normal  ovarian  surface  epithelial  cells  (IOSE)  or  epithelial  ovarian 
cancer  cells  (EOC)  using  a  Transwell  membrane  [41  ].  COLIIAI  mRNA 
was  statistically  significantly  upregulated  when  MSCs  and  INOFs 
were  co-cultured  with  EOC  but  not  IOSE  (Fig.  ID),  indicating  that 
cancer  cells  have  a  greater  capacity  than  normal  cells  to  induce 
COLIIAI  expression  in  fibroblasts. 

COL11A1  is  associated  with  cancer  progression  and  poor  survival 

COLIIAI  mRNA  expression  has  been  associated  with  poor  sur¬ 
vival  in  ovarian  cancer  [32,33]  and  kidney  cancer  [31  ].  To  elucidate 
the  underlying  biology  that  could  result  in  poor  survival,  we  inves¬ 
tigated  its  expression  in  ovarian  and  colon  cancers.  Using  a 
comprehensively  annotated  microarray  database  for  3431  human 
ovarian  cancers  [38],  we  show  that  increased  expression  of  COLIIAI 
mRNA  is  associated  with  overall  survival  and  disease-free  survival 
(Fig.  2A)  as  well  as  with  clinical  and  molecular  parameters  such  as 
increased  cancer  stage  and  grade  and  mesenchymal  molecular 
subtype  (Fig.  2B).  The  association  of  COLIIAI  expression  with  poor 
survival  is  unlikely  to  be  a  manifestation  of  the  total  amount  of 
stromal  fibroblasts  because  a  general  marker  of  fibroblasts,  vimentin 
{VIM),  is  not  associated  with  poor  survival  in  the  same  cohort  of 
ovarian  cancer  patients  (Table  SI ).  The  association  of  COLIIAt  with 


adverse  outcomes  is  also  not  restricted  to  ovarian  cancer.  We  show 
that  in  1820  colon  cancers  [42],  increased  expression  of  COLIIAI 
mRNA  is  associated  with  poor  disease-specific  and  disease-free  sur¬ 
vival  as  well  as  with  clinical  and  molecular  parameters,  such  as 
increased  cancer  stage  and  microsatellite  instability  and  CMS4  (mes¬ 
enchymal)  molecular  subtype  (Fig.  S1A,  B). 

To  systematically  investigate  an  association  between  COLIIAI 
mRNA  expression  and  survival  in  various  solid  and  liquid  cancers, 
we  plotted  COLIIAI  z-score  values  as  determined  by  the  pan¬ 
cancer  PREdiction  of  Clinical  Outcomes  from  Genomic  Profiles 
(PRECOG)  analysis  of -18000  cases  in  166  cancer  datasets  [36],  In 
most  epithelial  cancers,  COLIIAI  expression  was  associated  with 
poor  survival  (Fig.  3A).  Associations  between  expression  of  43  col¬ 
lagen  genes  and  survival  z-scores  in  12  common  epithelial  cancer 
types  revealed  that  for  the  majority  of  collagens,  increased  expres¬ 
sions  of  mRNA  were  associated  with  poor  survival,  with  COLIIAI 
having  the  strongest  association  (Fig.  3B). 

COL1 1  Al  is  among  the  most  differentially  expressed  genes  between 
cancers  and  corresponding  benign  tissues 

In  the  Genotype-Tissue  Expression  (GTEx)  project  database  [43], 
COLIIAI  mRNA  is  expressed  at  appreciable  levels  in  transformed 
skin  fibroblasts  but  not  in  non-transformed  skin  fibroblasts  or  other 
normal  tissues  (Fig.  S2).  Additional  analyses  of  various  expression 
datasets  containing  normal  adult  mouse  and  human  tissues 
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Fig.  2.  COLUAl  expression  is  associated  with  adverse  clinical  parameters.  (A)  Kaplan-Meier  analyses  of  overall  survival  (left  two  panels)  and  disease-free  survival  (right 
two  panels)  based  on  COLIIAI  expression  in  ovarian  carcinoma.  Disease-free  survival  includes  progression-  and  recurrence-free  survival.  Patients  were  stratified  to  COL11A1- 
high  (red)  or  COLllAl-low  (blue)  based  on  the  median  expression  of  COLUAl,  and  to  COL11A1-Q4  (highest  25%  expression;  red)  or  COL11A1-Q1  (lowest  25%  expression; 
blue).  (B)  COLUAl  expression  profiles  were  stratified  by  F1GO  stage  (left),  FIGO  grade  (middle),  and  molecular  subtype  (right).  Data  were  obtained  from  the  ovarian  microarray 
gene  expression  database  CSIOVDB  (csibio.nus.edu.sg/CSIOVDB/CSIOVDB.html).  HR,  hazard  ratio.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  this  article.) 
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Fig.  3.  COL11A1  expression  in  cancer  is  associated  with  poor  survival  in  multiple  cancer  types.  (A)  Survival  z-scores  in  different  cancer  types  associated  with  expression  of  COL11A1  mRNA.  (B)  Survival  z-scores  associated  with 
mRNA  expression  of  different  collagen  genes.  The  data  were  obtained  from  the  PREdiction  of  Clinical  Outcomes  from  Genomic  Profiles  (PRECOG)  database  (precog.stanford.edu). 
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revealed  that  COLiiAl  is  expressed  in  cartilage  and  collagen- 
producing  cells  in  the  eye  and  brain,  with  negligible  levels  in  most 
other  tissues  that  have  been  profiled  including  mesenchymal  stem 
cells  in  the  bone  marrow,  muscle,  and  fat  (data  not  shown). 

Comparison  of  COLiiAl  expression  in  17931  cancers  and  3503 
normal  tissues  (U133Plus2  platform)  and  9258  cancers  and 
4087  normal  tissues  (U133A  platform)  using  the  Gene  Expression 
across  Normal  and  Tumor  tissue  (GENT)  portal  [35]  revealed  that 
COLiiAl  mRNA  is  elevated  in  most  cancers  in  comparison  to  their 
corresponding  normal  tissues  (Fig.  4A).  In  some  cancers,  COLiiAl 
was  ranked  among  the  most  statistically  significant  differentially 
expressed  genes  when  cancer  and  its  corresponding  normal  tissue 
were  compared.  For  example,  comparison  of  cancer  and  normal 
tissues  in  The  Cancer  Genome  Atlas  (TCGA)  datasets  for  colon  cancer 
and  invasive  breast  cancer  ranked  COLiiAl  as  the  first  and  third  most 
differentially  expressed  gene,  respectively  (Fig.  S3). 

As  many  collagens  and  collagen-remodeling  genes  are  fre¬ 
quently  upregulated  in  fibroblast  activation  associated  with 
inflammation  and  fibrosis  in  the  absence  of  cancer,  use  of  these  genes 
as  therapeutic  targets  in  cancer  could  be  problematic.  Analysis  of 
expression  profile  datasets  show  that  levels  of  COLiiAl  mRNA  in 
inflamed  colonic  tissue  from  inflammatory  bowel  disease  and  fi- 
brotic  lung  tissues  are  not  significantly  different  from  those  in 
corresponding  unaffected  colon  and  lung  and  that  COLiiAl  levels 
associated  with  colon  inflammation  and  lung  fibrosis  are  minimal 
and  markedly  different  from  those  associated  with  colon  and  lung 
cancers  (Fig.  4B).  In  contrast,  levels  of  ACTA2,  the  gene  encoding  the 
prototypical  marker  of  myofibroblast  differentiation,  aSMA,  is  ex¬ 
pressed  at  similar  levels  in  cancers  and  inflamed  or  fibrotic  tissues 
(Fig.  4B). 

A  consistent  set  of  genes  is  co-expressed  with  COL1 1  Al  across 
different  cancers 

To  better  understand  the  biology  of  cancers  with  high  levels  of 
COLiiAl,  we  identified  genes  that  most  closely  correlate  with 
COLiiAl  mRNA  expression  in  13  TCGA  datasets  representing  dif¬ 
ferent  cancer  types.  Spearman’s  rank  correlations  between  COLiiAl 
and  its  co-expressed  genes  for  each  cancer  type  were  calculated. 
The  genes  were  then  ranked  based  on  the  average  correlation  of  each 
gene  across  the  13  cancer  types.  The  top  195  correlated  genes  were 
selected  based  on  an  average  correlation  of  >0.4  COLiiAl- correlated 
genes  were  then  ranked  based  on  the  average  of  the  absolute  cor¬ 
relation  values  (Table  1  and  Table  S2).  The  top  10%  most  highly 
correlated  genes  in  each  cancer  type  are  highlighted  in  pink  (Table  1 ). 
Notably,  COL  11 A 1  -correlated  genes  with  high  average  correlation 
scores  also  tended  to  be  among  the  top  10%  highest  scored  genes 
in  each  cancer  type  (indicated  in  pink  in  Table  1).  In  contrast,  the 
top  10%  COLllAl-anticorrelated  genes  were  not  conserved  across 
these  cancer  types  (Table  S3).  Some  of  the  top  ranked  COLllAl- 
anticorrelated  genes  in  individual  cancer  types  were  associated  with 
normal  functions  of  these  organs  suggesting  that  they  may  repre¬ 
sent  normal  tissue  or  a  noninvasive  tumor  component.  For  example, 
the  ovarian  cancer  top  100  COLllAl-anticorrelated  genes  (Table  S3) 
present  in  the  GSE121 72  ovarian  cancer  dataset  were  primarily  ex¬ 
pressed  in  ovarian  tumors  of  low  malignant  potential  (Fig.  S4). 

Pan-cancer  COL1 1A1  -correlated  genes  are  induced  in  CAFs 

Consistent  with  the  induced  expression  of  COL11A1  in  the  in  vitro 
co-culture  model  (Fig.  ID),  the  average  expression  of  the  pan¬ 
cancer  COL11A1- correlated  gene  set  was  significantly  induced  in 
mesenchymal  stem  cells  and  normal  ovarian  fibroblasts  co-cultured 
with  ovarian  cancer  cells  but  not  with  normal  ovarian  epithelial  cells 
(Fig.  S5).  Since  epithelial  cells  were  not  profiled  in  this  experi¬ 
ment,  it  is  unknown  if  fibroblasts  also  induce  expression  of  the  pan¬ 


cancer  COL  31A  3 -correlated  genes  in  epithelial  cells.  This  is  relevant 
because  several  of  the  195  pan-cancer  COLIIAI-c orrelated  genes  have 
been  shown  to  play  a  role  in  EMT  [42]  and  malignant  cells  under¬ 
going  EMT  have  been  proposed  as  one  possible  source  of  CAFs  [44], 
To  determine  if  the  pan-cancer  COL  3  3A  3 -correlated  gene  set  is  pref¬ 
erentially  expressed  in  cancer  cells  undergoing  EMT  or  in  host- 
derived  fibroblasts,  we  used  the  e-mtab-991  public  transcription 
profile  dataset  of  primary  patient-derived  colon  cancers  and  their 
patient-derived  xenografts  (PDX)  in  nude  mice  [39],  Presumably,  in 
PDX  samples,  fast-proliferating  human  cancer  cells  continued  to  grow 
in  mice  while  slow-proliferating  human  CAFs  were  lost  and  even¬ 
tually  replaced  by  mouse  fibroblasts,  which  can  be  distinguished 
from  human  cells  by  species-specific  gene  probes  [39].  GeneSet  clus¬ 
tering  analysis  showed  that  most  of  the  pan-cancer  COL11A1- 
correlated  genes  had  diminished  levels  in  PDX  samples  in  comparison 
to  primary  cancers  (Fig.  5A),  suggesting  that  the  genes  are  en¬ 
riched  in  the  CAFs  rather  than  in  the  cancer  cells.  However,  it  is  also 
possible  that  the  pan-cancer  COL  3  3A  3 -correlated  genes  are  ex¬ 
pressed  in  epithelial  cells  in  primary  colon  tumors  but  become 
silenced  upon  adaptation  of  human  cancer  cells  to  the  mouse  mi¬ 
croenvironment.  Thus,  we  conducted  GeneSet  clustering  analysis 
of  the  primary  ovarian  cancer  dataset  GSE40595  in  which  ovarian 
CAFs  and  epithelial  cancer  cells  were  isolated  by  laser-capture  mi¬ 
crodissection  [40].  The  pan-cancer  COLIIAI-c  orrelated  genes  were 
preferentially  expressed  in  CAFs  in  this  dataset  (Fig.  5B). 

In  addition  to  CAFs,  immune  cells  are  a  major  component  of  the 
tissue  microenvironment.  To  exclude  the  possibility  that  the  pan¬ 
cancer  COL33A3-correlated  gene  set  represents  immune  cells  in  the 
tumor  microenvironment,  we  used  the  expression  profile  of  230 
mouse  hematopoietic  cell  types  generated  by  the  Immunological 
Genome  Project  (ImmGen)  compendium  [45],  In  addition  to  he¬ 
matopoietic  cell  lineages,  the  dataset  contains  expression  profiles 
of  skin  fibroblasts  and  fibroblasts  residing  in  the  thymus,  lymph 
nodes,  and  spleen.  The  pan-cancer  COL1 1  Al  -correlated  gene  set  was 
highly  represented  in  fibroblasts  but  not  in  hematopoietic  cell  lin¬ 
eages  (Fig.  5C). 

CAFs  have  a  different  expression  profile  than  normal  fibro¬ 
blasts.  Moffitt  and  colleagues  defined  a  23-gene  signature  of ‘normal 
stroma’  and  a  25-gene  signature  of  ‘activated  stroma’  [46]  using  non¬ 
negative  matrix  factorization  for  virtual  microdissection  of  primary 
and  metastatic  pancreatic  ductal  cancer  samples  into  cell  subsets 
with  prognostic  and  biologic  relevance.  None  of  the  23  (0%)  ‘normal 
stroma’  genes  in  contrast  to  18  of  25  (72%)  ‘activated  stroma’  genes 
were  present  in  the  COL33A3-correlated  gene  set,  respectively 
(Fig.  5D),  suggesting  that  the  COL33A3-correlated  gene  set  repre¬ 
sents  CAFs. 

Biological  processes  associated  with  fibroblast  activation  in  cancer 

The  remarkable  uniformity  of  COLiiAl- correlated  genes  across 
13  different  cancer  types  suggests  involvement  of  these  genes  in 
common  biological  processes  that  are  independent  of  the  organ  site 
and  of  the  phenotypic  and  genetic  diversity  observed  in  individu¬ 
al  cancer  types.  To  gain  insight  into  the  biology  of  this  conserved 
gene  set,  we  conducted  several  analyses  that  identified  overlap 
between  the  195  pan-cancer  COL  3  3A  3 -correlated  genes  and  genes 
in  various  datasets  with  characterized  biological  features.  The  Gene 
Ontology  (GO)  Biological  Process  (BP)  analysis  revealed  that  the  pan¬ 
cancer  COLllAl-c orrelated  genes  are  primarily  involved  in 
extracellular  matrix  modification  and  collagen  remodeling  (Table  S4). 
Additionally,  we  used  SABiosciences  array  gene  tables,  which  consist 
of  literature-based  curated  molecular  pathways  where  each  pathway 
was  represented  by  84  genes.  Analysis  of  gene  overlap  between 
the  pan-cancer  COL  3  3A  3 -correlated  gene  set  and  genes  represen¬ 
tative  of  67  different  pathways  in  SABiosciences  arrays  showed  the 
largest  overlap  for  pathways  associated  with  extracellular  matrix, 
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Fig.  4.  Increased  expression  of  C0L11A1  in  cancer  and  low  expression  in  normal  tissues,  inflammation  and  fibrosis.  (A)  Comparison  of  C0L11A1  mRNA  expression  in  normal  tissues  and  corresponding  cancers.  Box  plots  for  two 
different  platforms  (U133Plus2  and  U133A)  were  generated  using  datasets  and  software  available  through  the  Gene  Expression  across  Normal  and  Tumor  tissue  (GENT)  portal  (medical-genome. kribb.re.kr/GENT).  The  y  axis 
shows  log2  mRNA  levels.  Average  expression  levels  in  normal  tissues  and  cancer  tissues  are  indicated  by  vertical  dotted  green  and  red  lines,  respectively.  (B)  C0L11A1  and  ACTA2  mRNA  expression  in  normal,  inflammatory  and 
fibrotic  conditions  in  comparison  to  cancer.  The  graphs  were  generated  using  the  public  R2  MegaSampler  software  (hgserverl.amc.nl/cgi-bin/r2/main.cgi)  for  the  processing  and  normalization  of  individual  datasets  imported 
from  the  Gene  Expression  Omnibus  (ul33p2,  MAS5.0  platform).  The  number  of  samples  in  each  GSE  dataset  is  indicated  in  parentheses.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred 
to  the  web  version  of  this  article.) 
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Table  1 

COL1 1A1  -correlated  genes  (Spearman’s  rank  correlation)  across  13  different  TCGA  carcinoma  types,  each  represented  by  >100  primary  tumors  from  therapy-naive  patients. 
Pink  denotes  the  top  10%  COL1 1A1  -correlated  genes  in  each  individual  carcinoma  type.  Rectangles  denote  genes  frequently  used  as  markers  of  activated  CAFs. 
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Fig.  5.  The  pan-cancer  COL  11  Al  -correlated  gene  set  is  expressed  in  CAFs  and  associated  with  poor  patient  survival  in  multiple  cancer  types.  (A)  GeneSet  expression  clustering  analysis  of  40  primary  colon  cancer  samples  and 
75  patient-derived  xenograft  (PDX)  samples  in  the  e-mtab-991  dataset  using  the  pan-cancer  COL11A1  -correlated  genes.  (B)  GeneSet  expression  clustering  analysis  of  laser-microdissected  ovarian  cancer  epithelial  cells  (32  samples) 
and  CAFs  (31  samples)  in  high  grade  serous  ovarian  cancer  in  the  GSE40595  dataset  using  the  pan-cancer  COL 11A 1  -correlated  genes.  The  Euclidean  distance  clustering  analysis  heatmaps  in  (A)  and  (B)  were  generated  using  the 
public  R2  GeneSet  Clustering  Analysis  tool  (hgserverl.amc.nl/cgi-bin/r2/main.cgi).  (C)  Expression  of  the  pan-cancer  COL  UA  1  -correlated  genes  mapped  on  the  transcriptome  of  individual  murine  hematopoietic  and  stromal  cell 
types  in  the  ImmGene  project  (immgen.com).  The  plot  was  generated  using  MyGeneset  tool  (rstats.immgen.org/MyGeneSet).  Black  dots  represent  mRNA  levels  (y  axis)  of  186  pan-cancer  C0L11A1  -correlated  genes  (9  genes 
were  not  present  in  the  database)  across  230  individual  cell  types  (X  axis)  grouped  into  10  main  groups.  (D)  Overlap  of  the  195  pan-cancer  C0LUA1  -correlated  genes  with  23  ‘normal  stroma’  and  25  ‘activated  stroma’  genes 
defined  by  Moffitt  et  al.  (E)  Unweighted  meta  z-scores  of  191  COL11A1  -correlated  genes  (4  genes  were  not  available  in  the  PRECOG  database)  were  compared  with  those  of  all  genes  (total  23,287  genes)  in  the  PRECOG  database 
using  unpaired  t  test.  The  plot  was  generated  using  GraphPad  Prism  software  version  6.0.  Intergroup  differences  were  assessed  by  the  Student’s  t-test.  Mean±SEM  of  pan-cancer  COL11A1  -correlated  genes  (0.8916  ±0.1641 
N  =  191);  mean  ±  SEM  of  all  genes  (0.09918  ±  0.01416  N  =  23287);  ****p<  0.0001. 
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fibrosis,  osteogenesis,  wound  healing,  EMT,  cardiovascular  disease, 
and  transforming  growth  factor  (3  (TGFfS)  signaling  (Table  S5). 
COL1I/U -correlated  gene  set  enrichment  analysis  of  chemical  and 
genetic  perturbations  (CGP)  showed  the  most  significant  overlap  with 
genes  up-regulated  in  association  with  cancer  invasiveness,  ad¬ 
vanced  stage,  stromal  cell  sternness,  and  epithelial-mesenchymal 
transition  (EMT)  (Table  S6).  The  uniformity  of  the  COL)  M2 -correlated 
genes  across  different  cancers  might  also  indicate  that  these  genes 
are  regulated  by  a  common  mechanism.  Ingenuity  Pathway  Anal¬ 
ysis  showed  that  transforming  growth  factor  beta  1  (TGFB1 )  is  the 
most  strongly  associated  upstream  regulator  of  the  pan-cancer 
COLU A) -correlated  genes  (Table  S7). 

COL11A1  -correlated  genes  are  associated  with  poor  patient  survival 
and  represent  potential  therapeutic  targets 

To  determine  whether  the  pan-cancer  COL  11 A I -correlated  gene 
set  is  associated  with  patient  survival  in  the  -18,000  cases  of  liquid 
and  solid  malignancies  in  the  PRECOG  dataset  [36],  we  compared 
survival  z-scores  for  the  195  pan-cancer  COL  11  Al -correlated  genes 
with  the  survival  z-scores  for  all  genes  in  the  dataset.  This  analy¬ 
sis  showed  that  expression  of  the  pan-cancer  COLllAl-correlated 
gene  set  is  significantly  associated  with  poor  survival  (Fig.  5E). 

Expression  profile  analyses  have  identified  a  mesenchymal  mo¬ 
lecular  subtype  of  cancer  associated  with  poor  survival  in  multiple 
cancers  including  ovarian  [47],  pancreatic  [46],  gastric  [48]  and  colon 
[49j.  In  colon  cancer,  it  has  been  shown  that  the  mesenchymal  mo¬ 
lecular  subtype,  which  constitutes  approximately  23%  of  colon 
cancers,  has  no  significant  enrichment  for  targetable  mutations  or 
copy  number  changes  in  candidate  driver  genes  [49  j.  Even  if  future 
research  identifies  targetable  events  in  cancer  cells  of  the  mesen¬ 
chymal  subtype,  it  is  predicted  that  enrichment  in  CAFs  and  excessive 
ECM  deposition  will  reduce  therapeutic  efficacy  by  creating  a  phys¬ 
ical  barrier  for  drug  transport.  Thus,  simultaneous  targeting  of  CAFs 
and  cancer  cells  may  be  necessary  for  chemotherapeutic  accessibility. 

To  identify  therapies  that  preferentially  target  activated  CAFs 
and  spare  normal  tissues,  we  combined  drug  target  searches  with 
expression  profile  datasets  in  cancers  and  normal  tissues.  Ingenu¬ 
ity  Pathway  Analysis  and  searches  of  the  Clinicaltrials.gov 
(clinicaltrials.gov)  and  ChEMBL  [50]  (ebi.ac.uk/chembl)  databases 
revealed  that,  of  the  195  pan-cancer  COL11A1 -correlated  genes,  16 
are  targets  of  drugs  used  in  clinical  trials  (Table  S8)  and  30  are  targets 
of  bioactive  compounds  (Table  S9). 

To  test  whether  any  of  the  drug/bioactive  compound  target  genes 
in  Tables  S8  and  S9  are  exclusively  expressed  in  activated  CAFs,  we 
determined  the  expression  levels  of  each  gene  in  normal  tissues  in 
the  GTEx  database  [43]  and  in  normal  vs.  cancer  tissues  in  the  GENT 
database  [35],  Additionally,  to  test  whether  selected  genes  were  ex¬ 
clusively  overexpressed  in  cancer  tissues  and  not  in  non-cancer 
associated  pathologies  such  as  inflammation  and  fibrosis,  we  com¬ 
pared  expression  of  these  genes  in  normal  tissues,  inflamed/ 
fibrotic  tissues  and  cancer  tissues  of  the  colon  and  lung.  Unlike 
COL11A1,  which  has  restricted  expression  in  normal  tissues  (Fig.  S2) 
and  is  highly  elevated  in  cancer  vs.  normal  tissues  (Fig.  4A)  and  in 
cancer  vs.  inflamed/fibrotic  tissues  (Fig.  4B)  most  of  the  target  genes 
were  expressed  at  high  levels  in  at  least  one  normal  tissue  and/or 
exhibited  equivalent  expression  levels  in  cancers  vs.  normal  tissues, 
and  cancers  vs.  inflamed/fibrotic  tissues.  One  example  of  this  pattern 
of  expression  is  CTCF  (Figs.  S6  and  S7).  However,  FN1,MMP13,  MMP14, 
FAP,  LOX  and  C0L1A2  exhibited  restricted  expression  in  normal  tissues 
and  elevated  expression  in  cancer  vs.  normal  tissues.  One  example 
of  this  pattern  of  expression  is  FAP  (Figs.  S8  and  S9).  Among  FN1, 
MMP13,  MMP14,  FAP,  LOX  and  COL1A2,  FAP  was  also  differentially 
expressed  between  inflamed/fibrotic  tissues  and  cancer  tissues  al¬ 
though  this  difference  in  expression  was  not  as  prominent  as  for 
COL11A1  (compare  Fig.  S9B  and  Fig.  4B). 


Discussion 

Whereas  in  the  past  most  therapeutic  approaches  have  focused 
on  the  cancer  cell  and  its  genetic  alterations,  it  is  becoming  appar¬ 
ent  that  the  microenvironment  plays  an  equally  important  role  in 
cancer  evolution.  We  now  recognize  that  the  cancer  stroma  not  only 
serves  as  a  scaffold  for  tissue  organization  and  integrity  but  also  pro¬ 
vides  key  biomechanical  and  molecular  signals  that  can  affect  various 
aspects  of  cancer  growth  and  biology,  including  proliferation,  sur¬ 
vival,  metabolism,  stem  cell  fate,  and  response  to  chemotherapy 
[51,52],  As  the  genetically  stable  subpopulations  of  the  cancer  mi¬ 
croenvironment  are  increasingly  recognized  as  potentially  effective 
therapeutic  targets,  a  comprehensive  definition  of  their  molecular 
characteristics  will  be  a  prerequisite  for  the  development  of  more 
precise  and  less  toxic  therapies.  Currently,  there  are  no  reliable 
methods  to  distinguish  activated  CAFs  from  non-activated  CAFs, 
which  although  frequently  abundant  within  cancers  do  not  neces¬ 
sarily  contribute  to  adverse  outcome.  We  identified  COL11A1  among 
the  top  differentially  expressed  genes  in  multiple  cancer  types  when 
cancer  tissues  and  their  corresponding  normal  tissues  were  com¬ 
pared.  We  showed  that  an  increase  in  COL11A1  expression  is 
associated  with  progression  and  poor  survival  in  most  cancer  types. 
COLUA1  is  a  particularly  attractive  therapeutic  target  because  of  its 
restricted  expression  in  normal  tissues  and  non-cancer  condi¬ 
tions,  such  as  inflammation  and  fibrosis. 

The  identification  of  a  highly  conserved  set  of  genes  associated 
with  COL11A1  expression  in  breast,  lung,  pancreas,  stomach,  urinary 
bladder,  colon,  thyroid,  cervix,  head  and  neck,  thyroid,  ovary,  and 
prostate  cancers  was  somewhat  surprising  in  light  of  the  genetic 
and  phenotypic  diversity  among  these  cancer  types.  The  con¬ 
served  expression  signature  indicates  that  the  reaction  of  stromal 
tissues  to  invading  epithelial  cancer  cells  may  be  similar  regard¬ 
less  of  the  organ  of  origin  or  genetic  alterations.  This  has  significant 
implications  for  the  development  of  pan-cancer  therapeutic  strat¬ 
egies.  Our  analysis  of  potential  upstream  regulators  of  the  pan¬ 
cancer  COL  11  Al -correlated  genes  revealed  TGFB1  as  the  most  likely 
candidate.  Dysregulation  of  TGF[3  signaling  is  recognized  as  the  main 
driver  of  fibroblast  activation  and  represents  the  most  logical  ther¬ 
apeutic  target  [53],  In  immortalized  normal  ovary  fibroblast  cell 
culture,  recombinant  TGFB1  has  been  shown  to  upregulate  expres¬ 
sion  of  COL11A1  and  several  other  COUJAl-correlated  genes;  this 
effect  was  abrogated  by  the  TGF|3  receptor  inhibitor  A83-01  [32]. 
However,  the  pleiotropic  nature  of  TGFfS  signaling  carries  the  risk 
of  adverse  effects  in  patients  [54].  In  order  to  abrogate  fibroblast 
activation  without  the  negative  effects  of  pan-TGFfS  therapy,  it  will 
be  necessary  to  design  therapies  for  more  specific  targets. 

Sixteen  of  the  195  pan-cancer  COLllAl-correlated  genes  are 
targets  of  drugs  in  clinical  trials.  These  targets  include  CTGF,  a 
matricellular  protein  involved  in  myofibroblast  formation  in  cancer 
as  a  binding  factor  of  fibronectin  and  a  downstream  mediator  of 
TGFfS.  A  clinical  trial  (clinicaltrials.gov;  NCT02210559)  is  currently 
enrolling  patients  with  unresectable  pancreatic  cancer  to  test  a  com¬ 
bination  of  conventional  chemotherapy  and  FG-3019,  the  human 
monoclonal  antibody  that  interferes  with  the  action  of  CTGF.  Our 
expression  analyses,  consistent  with  the  published  literature  (re¬ 
viewed  in  Ref.  [55]),  show  that  CTGF  is  expressed  at  similar  levels 
in  normal  tissues  and  cancers  and  therefore  unlikely  to  be  a  safe 
therapeutic  target  for  cancer  treatment.  In  contrast,  targets  such  as 
FN1,  FAP,  MMP13,  LOX  and  COL2A1  are  markedly  increased  in  cancer/ 
inflammation/fibrosis  compared  to  normal  tissues  and  are  thus 
predicted  to  have  a  better  safety  profile  than  agents  targeting  CTGF. 
Our  assessment  is  consistent  with  the  published  moderate  and  re¬ 
versible  toxicity  of  the  FN1  -targeting  monoclonal  antibody-cytokine 
fusion  protein  L19-IL2,  which  is  in  a  phase  I/II  study  for  patients 
with  solid  cancers  (clinicaltrials.gov;  NCT01058538)  [56,57],  Our  ex¬ 
pression  analyses  show  that  one  of  the  targets  of  bioactive 
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compounds,  FAP,  has  low  expression  in  normal  tissues  and  also  lower 
expression  in  inflamed/fibrotic  tissues  than  in  cancer.  Yet,  this  dif¬ 
ference  in  expression  may  not  be  sufficient  for  specific  targeting  of 
activated  CAFs  as  studies  in  mouse  models  have  shown  that  de¬ 
pletion  of  the  FAP+  stroma  can  induce  toxicity  due  to  expression 
of  FAP  in  the  mesenchymal  cells  of  bone  marrow,  muscle  and  adipose 
tissue  [58,59].  Future  efforts  to  specifically  target  activated  CAFs  can 
be  improved  by  designing  novel  therapies  to  target  genes  that  exhibit 
restricted  expression  in  nonmalignant  tissues.  When  considering 
COL11A1  as  a  cancer-specific  biomarker  and  therapeutic  target,  it 
is  important  to  note  that  several  normal  tissues  express  COL11A1. 
The  potential  side  effects  of  COL11A1  targeting  can  be  predicted 
based  on  the  phenotypes  of  mice  and  humans  expressing  mutant 
nonfunctional  forms  of  COL11A1.  A  homozygous  truncating  muta¬ 
tion  of  COLUA1  in  mice  results  in  poorly  formed  cartilage  [60], 
while  human  COLltAl  mutations  are  associated  with  articular 
hypermobility,  dermal  hyperelasticity  and  widespread  tissue  fra¬ 
gility  [61].  Of  note,  these  collagenopathies  are  associated  with  the 
absence  of  COL11A1  function  throughout  development  and  are  un¬ 
likely  to  manifest  upon  transient  targeting  of  COL11A1  in  adults. 
Additionally,  since  COL11A1  and  many  of  the  pan-cancer  COL11A1- 
coexpressed  genes  have  multiple  tissue-specific  mRNA  splicing 
isoforms,  it  will  be  valuable  for  future  targeting  purposes  to  deter¬ 
mine  if  any  mRNA  isoforms  are  specifically  expressed  in  activated 
CAFs  [62], 
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